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Abstract 

Fractional calculus is the calculus of differentiation and integration of non-integer 
orders. In a recent paper (Annals of Physics 323 (2008) 2756-2778), the Funda- 
mental Theorem of Fractional Calculus is highlighted. Based on this theorem, in 
this paper we introduce fractional series expansion method to fractional calculus. 
We define a kind of fractional Taylor series of infinitely fractionally-differentiable 
functions. By using this kind of fractional Taylor series, we give a fractional general- 
ization of hyper geometric functions and derive corresponding differential equations. 
For finitely fractionally-differentiable functions, we observe that the non- infinitely 
fractionally-differentiability is due to more than one fractional indices. We expand 
functions with two fractional indices and illustrate how this kind of series expansion 
can help to solve fractional differential equations. 
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1 Introduction 



Fractional calculus is the calculus of differentiation and integration of non-integer orders. 
During last three decades or so, fractional calculus has gained much attention due to its 
demonstrated applications in various fields of science and engineering [IJ[2j[3l[Il[5l[6j[7J[El 
El HOI HH H2l [131 [HI H5l UHl H71 H8] . There are many good textbooks of fractional calculus 
and fractional differential equations, such as [H El El IU E]- For various applications of 
fractional calculus in physics, see [SI El El 13 HDl HH H21 H31 H31 H51 HSl IIZl HB| axid. references 
therein. 

In Calculus, the Fundamental Theorem of Calculus (Newton-Leibniz Theorem) is of 
fundamental importance, matching its name. For fractional calculus, an analogous theo- 
rem has been highlighted recently in a paper |18j . When trying to construct a consistent 
fractional vector calculus, V. E. Tarasov observed that many of fundamental problems 
can be solved by using a generalization of the Fundamental Theorem of Calculus. 

Series expansion is an important tool to calculus. Particularly, series expansion plays 
an important role in solving some differential equations, such as the hypergeometric dif- 
ferential equations [T9l [20l |2~T] . However, fractional series expansion has not yet been 
introduced to fractional calculus. 

In this paper, based on the Fundamental Theorem of Fractional Calculus, we will in- 
troduce fractional series expansion method to fractional calculus. We will define a kind 
of fractional Taylor series of infinitely fractionally-differentiable functions. By using of 
this kind of fractional Taylor series, we will give a fractional generalization of hypergeo- 
metric functions and derive corresponding differential equations. For finitely fractionally- 
differentiable functions, we observe that the non-infinitely fractionally-differentiability is 
due to more than one fractional indices. We will expand functions with two fractional 
indices and illustrate how this kind of series expansion can help to solve fractional differ- 
ential equations. 

The structure of this paper is as follows. In Section El we briefly review fractional 
derivative, fractional integral and the Fundamental Theorem of Fractional Calculus. In 
SectionEl we introduce the fractional Taylor series of an infinitely fractionally-differentiable 
function and give some examples. In Section HI we make a generalization of the hyper- 
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geometric functions. In Section [51 we discuss finitely fractionally-differentiable functions. 
In Section^, we give our summary. 



2 Fractional Caculus 

In this section, we briefly review the definitions of fractional integral and fractional 
derivative, and the Fundamental Theorem of Fractional Calculus. For more details, see 

[H [21 El HE]. 

2.1 Fractional integral and fractional derivative 

There are many ways to define fractional derivative and fractional integral. Most of them 
are based on the idea that we can generalize the equation 

d n x m ml 



dx n (m — n) 
by replacing each factorial by a gamma function, to 

fi" r(p + i 



x m ~ n , neN, 1 



/3 " a , a > 0. (2) 



dx a T(f3-a + 1) 
In terms of integral operation /, the idea is to generalize 

777 1 

Fx m = —x m+n (3) 

(m + n)l v ' 

to 

Fx? = w r(/3 t 1} y +a . (4) 

The most commonly-used fractional integral is the Riemann-Liouville fractional inte- 
gral (RLFI). 

Let f(x) be a function defined on the interval [a, b]. Let a be a positive real. 
The right Riemann-Liouville fractional integral (right-RLFI) is defined by 

Ia\J{x) = p^y f {x - O^/CO*, (5) 
and the left Riemann-Liouville fractional integral (left-RLFI) is defined by 

i: lb m = j m - xr-'fim- (6) 



We note that in this paper our use of notations "right" and "left" is different from the 
common use, for reasons that will be clear later. 

For fractional derivatives, the Caputo fractional derivative (CFD) is a commonly-used 
one. Let n = [a] + 1. 

The right-CFD and left-CFD are defined, respectively by 

One can check that the above definitions really generalize ([1]) and ([3]), and give 

D C ^(x-af = r( ^^ ) 1) (x-a/-°, + 0, 1, .... [a]; (9) 
D°g{b-x)l> = r^-a+l) ^"^""' ^O' 1 ' -.M; < 10 > 

= r^TiJ*-^- (12) 

Especially, the Caputo fractional derivative on a constant (/3 = 0) yields zero, 

D C a \ x a -1 = 0, (13) 

D C 4 -1 = 0. (14) 

This simple property is decisive in the fractional series expansion and in our preference 
of the Caputo fractional derivative to another commonly-used fractional derivative, the 
Riemann-Liouville fractional derivative, whose operation on a constant gives not zero, 

D RL > a -l= ( x ~ a ^ a <a <l, (15) 
«l* r(l-a)' V ' 

D ^' a ■ 1 = < a < 1. (16) 
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2.2 Fundamental Theorem of Fractional Calculus 

The Fundamental Theorem of Calculus (Newton-Leibniz Theorem) is 

J b dx ±f( x ) = f(b)-f(a), (17) 

f IiO<K /(*). (18) 

J a 

This theorem means that the derivative operation is inverse to the integral operation, and 
vice versa. 

In fractional calculus, an analogous theorem exists [51 [T8]. 
Fundamental Theorem of Fractional Calculus. 

(1) Let a > and let f(x) G L^a, b) or f(x) G C[a, b]. Then 

D^%f(x) = f(x), (19) 
D C x p:\ b f{x) = /(*). (20) 

(2) Let < a < 1. If /(x) G AC[o, 6] or f(x) G C[a, 6]. Then 

i: ]x D^f(x) = f(x)-f(a), (21) 
= f(x)-f{b). (22) 

Here Loo(a, 6) is the set of those Lebesgue measurable functions / on (a, b) for which 
H/lloo < °° j where H/Hoo = ess sup a < :r <ft|/(a;)|. Here ess sup|/(x)| is the essential max- 
imum of the function AC [a, b] is the space of functions / which are absolutely 
continuous on [a, b]. C[a, b] is the space of continuous functions / on [a, b] with the norm 
\\f\\c = max x . eM 

The proof of this theorem can be obtained from [5], in which these results are included 
in Lemma 2.21 and Lemma 2.22 there. 

So, by this theorem, one can say that the right (left) Caputo fractional derivative 
operation and the right (left) Riemann-Liouville fractional integral operation are inverse 
to each other. 
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3 Fractional Taylor series of infinitely fractionally- 
differentiable functions 



In this section, we will introduce fractional series expansion method to fractional calculus 
and define a kind of fractional Taylor series. 
We observe that 

/(*) = f(a) + D a al J(x)\ x= ^ ■ (I% x • 1), (23) 

where a < £ < x, and £ varies with the integral upper bound. Here and after, D^ x denotes 
the right Caputo fractional derivative D,'°. 

The corresponding formula in integer-order calculus is 



/<*)=/<-)+! 

which is the Lagrange Mean Value Theorem. 
Make a step further, 



■(x-a), (24) 



f(x) = f(a) + D a al J(x)\ x=a ■ (i: ]x • 1) + D^D^f(x)\ ■ (1^ ■ 1). (25) 



The corresponding formula in integer-order calculus is 



•(-«) + (£) 2 / 

x=a 



■\{x-af. (26) 



And so on. One can extend this expansion to infinite order if the function is sufficiently 
smooth. 

Based on this observation a definition of a formal fractional Taylor series expansion 
can be made. 

Definition l.a. Let f(x) be a function defined on the right neighborhood of a, and be 
an infinitely fractionally-differentiable function at a, that is to say, all (D^ x ) m f(x)(m = 
0, 1,2,3..) exist, and are not singular at a. The formal fractional right-RL Taylor series 
of a function is 

oo 

/(*) = J2( D «\*) m fWl=a ■ ■ !]■ (27) 

m=0 

Explicitly, 

(/« ) m .l = — J— ( x - a ) ma . (28) 
v alxJ T(ma + 1) 
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Definition l.b. Let f(x) be a function defined on the left neighborhood of b, and be 
an infinitely fractionally-differentiable function at b, that is to say, all (D^ b ) m f(x)(m = 
0, 1,2,3..) exist, and are not singular at b. The formal fractional left-RL Taylor series of 
a function is 

oo 

/(*) = E(^|6) m /(z)U • t(4V m ■ !]• ( 29 ) 

m=0 

Explicitly, 

(4V m -i = r , 1 -—A b-x) ma . (30) 

1 (ma + 1) 

In the above definitions, D^ x is the right Caputo fractional derivative D^£; D^ b is the 
left Caputo fractional derivative D^ b . I^ x and I£j 6 are right- and left- Riemann-Liouvelle 
fractional integral, respectively. 

Remark 1. One can easily check the formal correctness of the expansions by using of 
the Fundamental Theorem of Fractional Calculus, or the relations (l9])- ffT4"j) . For rigorous 
validity, convergence is required. 

Remark 2. Series expansion has played an important role in Calculus, particularly 
in solving differential equations. However, fractional series expansion has not yet been 
introduced to fractional calculus. This is because a pre-requisite that makes fractional 
series expansion possible is the Fundamental Theorem of Fractional Calculus, which is 
only recently proved and highlighted [5J [18] . 

Remark 3. We may expect that fractional series expansion will shed new light on 
fractional calculus, especially the field of fractional differential equations. In the next 
section, we will use this expansion to define a fractional generalization of hypergeometric 
functions and discuss their differential equations. 

Remark 4. The fractional Taylor series is defined for infinitely fractionally-differentiable 
functions. Finitely fractionally-differentiable functions will be discussed in Section |5j 

In the following, we give some simple examples of fractional Taylor series. 
Example 1. (x — a) 13 , with no / G iV satisfying (3 = la. 

W\T(x-«Y= T ^^ ] +l) (x-°Y- (3D 

For large m, the derivative will be singular at a. So we cannot make the expansion. 
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Example 2. 

m=0 ro=0 v y 

Example 3. The Mittag-Leffler function E a ((x — a) a ), which satisfies 

D% x E a ((x - a) a ) = E a ((x - a) a ), E a (0) = 1, (33) 

It is the fractional analogue of exp(x). For arbitrary to, (D^ x ) m E a ((x — a) a )\ x=a = 1. So, 

OO OO 

Ea ((x - a )°) = £[(ig x r • i] = E TwTT) (a; - a)ma - (34) 

m=0 m=0 ^ ' 

Notice the difference between e^~ a ^ and E a ((x — a) a ). 
Example 4. cos a ((x — a) a ) and sin a ((x — a) a ). 

cos (ix an 1 iX ~ ar ' 2 I (X ~ ar4 (X ~ ar I (35) 
cos a {{x a))-l r(a . 2 + 1) + r(a . 4 + 1) r(a-6 + l) + - {6b) 



(x — a) a ' 3 (x — a) a ' 



sin a ((x - a)*) = 1 ' - \ ' + v ' - ... (36) 
1 (a + 1) 1 (a ■ 6 + 1) 1 (a ■ 5 + 1) 



They satisfy 



i^sm Q ((a;-a) Q ) = cos a ((x - a) a ), (37) 
D% x cos a ((x - a) a ) = -sin a {{x-a) a ). (38) 

4 Fractional hypergeometric function 

In the section, we will define a fractional generalization of the hypergeometric functions. 

Let us first consider a fractional generalization of the confluent hypergeometric differ- 
ential equation: 

z a {D% g ) 2 y + (c - z a )D% ]z y -ay = 0. (39) 

Here a and c are complex parameters. When a — 1, this is the ordinary confluent 
hypergeometric equation. 

Introducing the fractional Taylor series 

OO 

y{z) = Y,c k z a -\ (40) 
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and substituting, we get the ratio of successive coefficients 

rfte+i) 



c k+1 ■ T(ka + a + 1) a + 



r(fca-Q!+l) 



c k -T(ka + l) r (te+ i) 

F(ka— a+1) 

c\ ■ Y{a. + 1) a 
c c' 

Thus we get a solution of the above differential equation, 



(41) 
(42) 



y {z) = Y s*°" k - (43) 



Here (a)£ is defined as 



[o)o = 1, («)? 



/ iYfcot - a + 1) \ / r(fca-2a + l)\ . ,„ 7 

<°K = ( a+ r ( L-2 a + 1 ) )( a+ r( to -3 a + i) )- (a>? - i£2 ' (44) 

This can be seen as a fractional generalization of the rising factorial 

(a) k = (a + k - l)(a + fc - 2). ..a. (45) 



And the series (J43J) can be seen as a generalization the confluent hypergeometric function. 
If a = 1, it is exactly the confluent hypergeometric function. 

For the fractional Gauss hypergeometric function, consider the following series 



{z) = y (a)g(6)g x k 

U (c)g T(A;a + 1) 



which reduces to the Gauss hypergeometric function when a — 1. 
The ratio of successive coefficients is 



c fc+i " T(/ca + a + 1) 
Cfc • r(fca + I) 



r(fca+l) \ / , , r(fca+l) 
u "r" r(to-a+l) / I ~ l ~ r(fca-a+l) 



c i r(fca+l) 



r(fca-a+l) 



(47) 



Making some manipulation, one can get 

V(ka + l) T(ka + a + l 



T(ka — a + 1) Y{ka — a + 1) 

h-L r x«JfclLx r(fca + i) r(fca+i) 

c fc -a6 + c fc -(a + 6)— — + c k ■ — -r^ — r- (48 

1 {ka — a + 1 [ka — a + 1) 1 [ka — a + 1) 



9 



This equation can be translated to a fractional differential equation 

ah ■ f(z) + (a + b)(z - z ) a D« M J{z) + (z - Zo) a D^ {M [(z - z^D^Jiz)] 

= c ■ D« M J{z) + (z- z r(Dl ]z ) 2 f(z). (49) 

When a = 1, this equation reduces to the ordinary Gauss hypergeo metric equation. One 
can check y(z — z ) defined in (|4T)|) satisfies this equation. 
For generalized fractional hypergeometric series 

making repeated use of (z — zo) a D^ z , one can also get its differential equation. 

5 Functions with two fractional indices 

Not all functions are infinitely fractionally-differentiable, so it is meaningful to investigate 
finitely fractionally-differentiable functions. In Example 1 in Section [31 we have given an 
example function that is not infinitely fractionally-differentiable. 

We observe that the non-infinitely fractionally-differentiability is due to another frac- 
tional index j3 (no / satisfying (3 = la). A function f(x) is said to have the behavior of 
fractional index a in the right neighborhood of a, if it can be expanded as a series with 
the basis {(x - a) ma \m = 0, 1, 2, 3...}. 

We also observe that some functions can be expanded into a form with two fractional 
indices, i.e. 

oo oo 
m=0 n=0 

These functions can reduce to a function that is infinitely a-differentiable or infinitely 
/3-differentiable, but generally they are finitely fractionally-differentiable functions. 

When we are solving fractional differential equations, we should take care if [a] > 
and (3 G N, for the reason of the relation (J9]). However, for [a] = or (3 not an integer, 
the above series is fairly a good ansatz. 



10 



A function f(x) is said to be N • a fractionally-differentiable at a, if (D^ x ) m f (x)\ x=a 
is finite for m < N, but infinite for m — N + 1. Then it could be expanded as 

f(x) = c 00 + c 10 (x-a) a + +c 20 (x-a) 2a + ... + c N0 (x - a) Na 

+ c mn (x - a) ma+n ? + c mn (x-a) ma+n P, (52) 

Na<ma+nf3<(N+l)a ma+nf3> (N+l)a 

in which the first term in the second line should not be zero. 

Now let us see how this kind of series expansion could help to solve fractional differ- 
ential equations. 

Consider the following equation: 



(* " *) a M| J 2 f(x) - D%J(x) = x _ x + a ' ■ (53) 
The right hand side can be expanded as: 

- r.i" a a) = s> - -) k + s> - - <■ ( 54 ) 

k=l k=0 

It is of index a and index (3—1. It is 1 • a fractionally-differentiable. So f(x) is 2 • a 
fractionally-differentiable. We can write 

f(x) = c 00 + c w (x - a) a + c 20 (x - a) 2a 

+ + c n (x - a) a+1 

+ c 02 {x - a) 2 

+ J2 c mn (x-a) ma+n . (55) 

ma+n>3a,(m,n)^(0,2) 

Here we assume < a < 1. 
Or more manageably, 

f(x) = c 00 + c w (x - a) a 

CO 

\a+n 



+ E ° 0n ( X ~ ^ + E Cl "( ;r ~ 
ra=2 ra=l 

CO CO 

+ c ^ x - a ) 2Q+n + E E c -"( x - a ) ma+n - ( 56 ) 



n=0 m>3 n=0 
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Substituting fl54|) and fl56l) into (1531 . one gets 

Cmn = 0, ?T7, > 3; 

T(n + 2a + l) r(n + 2a + l) 
r(n + 1) T(n + a + 1) ' 



1/C2 



+ a + + a + 

!/cm = ^7 — r p / — Tx , n > 1; 

r(n — a + 1) r(n + 1) 

Cio = 0; 

c n = 0, n > 1; c 00 = /(a)- (57) 
Thus we solved the fractional differential equation ( 1531) . 

The above procedure can apply generally to fractional differential equations. For a 
generic fractional differential equation F[(x — a),y(x — a),D^ x ] = 0, we can solve it by 
the above procedure summarized as follows. 

1. Expand the (x — a) part of the equation and find out the indices of y(x — a), 

2. Find out the number N such that y{x — a) is iV • a fractionally-differentiable, 

3. Expand y(x — a), 

4. Substitute the expansion series into the equation, 

5. Find out the coefficients of the expansion series of y(x — a). 

We give another example in the following. Consider the fractional differential equation: 

D^y{x)+f{x)y(x) = g(x), (58) 

where f(x) and g(x) are given functions of index /3, 

fix) = f Q + h ■ (x - af + f 2 -{x- a) 2 ? + f 3 ■ (x - a)^ + (59) 

9(x) = 9o + gi-{x- af + g 2 -{x- a) w + g 3 ■ (x - a) w + (60) 
where /j and g^ (i = 0, 1, 2, 3, ...) are constants. 
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So y(x) have two indices a and (3. If g\ 7^ and a < (3 < 2a, g(x) is 1 ■ a fractionally- 
differentiable, then y(x) is 2 • a fractionally-differentiable. We can write 

y(x) = c 00 + ci (x - a) a + c 20 (x - a) 2a 

+ + c n (x-a) a+/3 

+ c 02 (x-a) 2/3 

+ c mn (x-a) ma+n P. (61) 

ma+n)3>3«, (m,ra)^(0,2) 

Or more manageably, 

= c 00 + ci (aj - a) a 

00 00 

n=2 n=l 

oo 00 

+ ^c 2n (x-a) 2a+ ^ + ^^c mn (x-a) ma+ ^. (62) 

n=0 m>3 n=0 

Substituting (|55j) . p?|) and (|B2]1 into ([58]), one gets 

coo = y(a), c 0n = 0, n > 1; 

T(n/3 + l) 
cm - f(aT^+l) ( ^" /nC00); 

r(q + n/3 + l) 

r(2a + + 1) 

r((m - l)a + n/3 + 1) 

"C(m— l)n 



r(m« + n/3 + 1) 
xm _ 2 T(2a + n(3+l) 



1 (ma + np + 1) 
Thus we solved the fractional differential equation (1581) . 

In this section we have discussed functions with two fractional indices, but the exten- 
sion to functions with more fractional indices will not be difficult. 



6 Summary 

In summary, in this paper we introduced fractional series expansion method to frac- 
tional calculus. We defined a kind of fractional Taylor series of infinitely fractionally- 
differentiable functions. Based on our definition we generalized hypergeometric functions 
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and derived their differential equations. For finitely fractionally-differentiable functions, 
we observed that the non-infinitely fractionally-differentiability is due to more than one 
fractional indices. We expanded functions with two fractional indices and illustrated how 
this kind of series expansion can help to solve fractional differential equations. 
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